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Abstract: Symanzik effective actions, conjectured to describe lattice artifacts, are de- 

^ . termined for a class of lattice regularizations of the non- linear O(A^) sigma model in two 

dimensions in the leading order of the 1/A^-expansion. The class of actions considered in- 



CN \ eludes also ones which do not have the usual classical limit and are not (so far) treatable in 

^■'^ ' the framework of ordinary perturbation theory. The effective actions obtained are shown 



to reproduce previously computed lattice artifacts of the step scaling functions defined in 



CO ' finite volume, giving further confidence in Symanzik's theory of lattice artifacts. 
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1 Introduction 

Determinations of physical quantities in lattice QCD require extrapolation of the data to 
the continuum limit. For this purpose an ansatz must be made on the nature of the lattice 
artifacts. Usually these are assumed to be predominantly power-like ~ o^ in the lattice 
spacing a, the power p depending on the particular lattice action used for the simulation. 
The latter ansatz is based on a conjecture of Symanzik [1] which says that the lat- 
tice artifacts of correlation functions of the basic fields are described by a local effective 
continuum Lagrangian involving higher dimensional operators which may have only the 
symmetries of the underlying lattice theory, and coefficients depending on a ^ . 



^The description of lattice artifacts of correlation functions involving composite operators involve in 
addition further effective operators. 



The study of lattice artifacts is not only of theoretical interest: it is closely related to 
Symanzik's improvement program [2-5] where the leading cutoff effects are cancelled by 
adding (already on the lattice) O (a^) local counterterms with suitably fine-tuned coeffi- 
cients. Using such improved actions the lattice artifacts are reduced considerably and since 
the continuum limit is reached faster, allows for cheaper MC simulations. 

This structural conjecture is based mainly on perturbative computations, both in the 
coupling and in the 1/A^ expansion, in various models. An all order proof for the exis- 
tence of Symanzik effective actions for 01 theory and QED4 has been given by Keller [6], 
albeit this proof is only in the framework of a continuum regularization. However, these 
theories, at least when regularized on the lattice, are probably trivial in the continuum 
limit. Renormalized couplings tend to zero as g ~ c/ ln(a/u) and hence the continuum limit 
(a free theory) is actually reached only logarithmically! Treating a renormalized coupling 
effectively as a constant for a range of cutoffs one has for small g a perturbative Lagrangian 
description of the low energy physics, and in this case the Symanzik effective Lagrangian 
describes the leading a} cutoff corrections to this, in particular effects of rotation symmetry 
breaking. 

So far an all order perturbative proof of the conjecture has not been given for asymp- 
totically free theories; there seems to be no technical obstacle but it is probably a tedious 
task. For QCD there is support for the validity of the conjecture coming only from low order 
perturbative computations, and from the fact that the numerical data can be convincingly 
described by fits of the expected form. 

Asymptotically free theories in 2 dimensions such as the O(A^) sigma model, offer 
a laboratory to study cutoff effects in more detail, since simulations to large correlation 
lengths can be made, and one can also perform studies in the 1/A'^-expansion. We shall 
restrict attention to this theory in this paper and the structure of the Symanzik effective 
action is discussed in the next subsection. For a class of lattice actions permitting a 
standard perturbative formulation it was shown [7] that the generic leading cutoff behavior 
is oc a^(ln?n,a) '' ^■^^ (see subsection 1.2). The additional In a factor for A^ = 3 becomes 
visible in precision measurements, and can even mimic O (a) effects in some intermediate 
range of correlation lengths. 

In a recent paper [8] a class of models was identified with very small lattice artifacts. 
Some of these models are particularly fascinating since they seem to belong to the same 
universality class as the standard-type actions described above but don't have the usual 
classical limit and so far a perturbative treatment for them is unknown (see subsection 1.2). 
However, they can be treated in the 1/A approximation and there the step scaling functions 
(of the finite volume mass gap) were found to be of the form c? In^ a with s = 2 in leading 
order 1/A rather than lim7v->oo N/{N — 2) = 1. Note that in the framework of Symazik's 
conjecture, the difference between the cutoff behavior for different lattice actions can only 
arise from the couplings in the effective theory. 

In this paper we determine the Symanzik effective action for the classes of lattice actions 
mentioned above in the leading order of the 1/A expansion by matching lattice correlation 
functions to the effective theory. We show how the constructed effective theory reproduces 
the artifacts of the step scaling functions defined in finite volume. It is instructive to see 



how the additional non-perturbative effects leading to the In a appear in this framework. 
The various steps of the computation are given in sections 2-4 and some technical details 
are relegated to appendices. 

1.1 Symanzik's effective action for 2d 0{N) sigma models 

As explained above, Symanzik's idea is to mimic the lattice artifacts (at least up to the 
leading O (a^) order) by using an effective Lagrangian in the continuum theory: 

-C,s = -C + a^^Ciig)Ui, (1.1) 

i 

where C is the original continuum Lagrangian density and in this sum all operators Ui 
with engineering dimension 4 and reflecting the lattice symmetries have to be taken into 
account. The c-number coefficients Ci{g) depend on the lattice action. For perturbative 
type actions they can be calculated in perturbation theory. In this paper we will calculate 
them in the large N expansion. 

For practical calculations a useful starting point is the master formula of Symanzik's 
effective theory, which can be obtained using the effective Lagrangian (1.1) and is written 
directly in terms of correlation functions on the lattice and the corresponding effective 
continuum model: 

g^,,,{Xo,a)=f{g)gf^^{g,a-') + a'f{g)Y,v^{9)gf^n)i9,a-') + 0{a^) . (1.2) 

i 

Here G^^^ represents a general lattice correlation function (with r external legs) for any 
physical quantity X, in real space or momentum space, and G^j^^ is the analogous renor- 
malized quantity calculated in the continuum model. G^^^ depends on the set Aq of bare 
lattice couplings and parameters, while the continuum correlators depend on the renormal- 
ized coupling g and the renormalization scale /i, which is taken here, for simplicity, as the 
inverse lattice spacing a^^. G^n\ is the corresponding renormalized continuum correlation 
function calculated with insertion of the integral of one of the local dimension 4 operators, 
Ui, appearing in Symanzik's effective continuum action. For a given lattice symmetry the 
set of operators is fixed and only the c-number coefficient functions Vi{g), the renormalized 
coefficients corresponding to Ci{g), depend on the lattice action. Finally, the wave function 
renormalization constant y{g) and the relation g = g{Xo) between the lattice coupling pa- 
rameters and the continuum coupling constant are again action dependent. The latter can 
be obtained by calculating some physical mass parameter Ai on both sides: 

A1(g,a-i)=Matt(Ao). (1.3) 

For generating correlation functions we will use the source dependent action 

A= /'d^x£= fd'^xl^d^S-d^S-^I-sl . (1.4) 

Here S is the A^-component sigma model field with normalization S-S = 1 and instead of the 
bare coupling constant g^ we will mainly use the 't Hooft coupling /o = Ng^ (and similarly 



for the renormalized quantities: / = Ng^). In the continuum we wih use dimensional 
regularization va. D = 2 — e dimensions. 

The set of local operators appearing in Symanzik's effective action always includes the 
following Lorentz-scalar dimension 4 O(A^) invariant operators: 



(1.5) 



These operators mix under renormalization with other Lorentz-scalar operators and in 
Symanzik's effective action we also have to include operators which are not fully Lorentz 
scalar but invariant under lattice symmetries only. The full list of operators in the case of 
lattices with cubic symmetry will be given in section 2. 

A set of lattice regularized 0{N) sigma models with one tunable coupling /? = 1/Aq is 
given by the following quadratic lattice actions: 
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^latt = f «' E '^(^) • ^(y)^(^ - y) ■ (1-6) 

x,y 

Here the only requirement is that the Fourier transform of the inverse "propagator" K(x) 
defined by 

Kp = a2^e-^"P^i^(x) (1.7) 

X 

behaves for a — t- as 

Kp=p'^ + (a^) . (1.8) 

This ensures that the lattice model has the correct classical continuum limit. For the 
standard lattice action (ST) we have: 

An alternative way of writing the master equation is 

eiatt = ^''^°^(Ao,a){l + a2<5''(Ao,a) + 0(a4)} , (1.10) 

where the scaling part Q'^'^^> is universal (up to wave function renormalization) and the 
correction factor is of the form 

5''iXo,a) = Y,v^i9)Sfi9,a). (1.11) 

i 

Using perturbation theory and renormalization group considerations, it was shown in [7] 
that close to the continuum limit this correction can be represented as 

5^(Ao,a) = CiDf(A){/3i+2x + o(/32x)} + o(;3>^) , (1.12) 



where x = l/(-^ ~ 2) and close to the continuum 

/3 = ^ ~ In(Aa) . (1.13) 

In the above formula Ci is a non-universal (action dependent) constant which is calculable 
in perturbation theory. On the other hand, Di{K) is universal but non-perturbative and 
depends on the quantity X in question and on the non-perturbative physical scale A. 
D^ {A.) is related to the non-perturbative matrix element of the operator with largest 
anomalous dimension. For A^ = 3 (1-12) predicts a large logarithmic correction proportional 
to a^(lna)^, whereas for large A^ this prediction gives O(a^lna) corrections. This is 
consistent with the known [9, 10] O (a^ In a) corrections appearing for the standard lattice 
action in the large A^ limit. 

1.2 Lattice artifacts for non-perturbative actions 

The universality class of lattice regularized 0{N) sigma models is actually much larger than 
the set of quadratic actions (1.6). In [8] non-perturbative actions were studied. For these 
types of actions the classical continuum limit cannot be expanded by usual weak coupling 
perturbation theory or there is no classical limit at all! Nevertheless, as was demonstrated 
in [8] the exact quantum continuum limit of these models coincides with the usual one. 
Moreover, some of these models show much smaller lattice artifacts than the conventional 
ones described by (1.6) type actions. 

The physical quantity studied in [8] was the step scaling function cr(2, u) which is 
defined [11] as follows. Consider the model confined in a periodic (one-dimensional) box 
of size L and denote by M{L) the mass gap in the box.^ Introducing the dimensionless 
LWW couphng [11] 

u = LM{L), (1.14) 

the step scaling function 

a{2,u)=2LM{2L) (1.15) 

describes how the LWW coupling changes under doubling the size of the box (a discrete 
renormalization group transformation). On the lattice we have 

u = LM{L,a), S(2,n,a/L) = 2LM(2L,a). (1.16) 

We will also use the notations 

u' = T.{2,u,a/L), u'^ = a{2,u) = J:{2,u,0). (1.17) 

For large A^ the step scaling function is known exactly [12]. It is given by the solution 
of the implicit equation 

/o(n) = /o(n:^) + ^ . (1.18) 

The function fo{u) will be given explicitly by (2.74). 



^M or M(L) always refers here the finite-volume mass gap, while we denote the infinite-volume mass 
gap by m. 



Let us introduce the variable^ 



i-l„[^V (1.19) 

2tt \ ma ^ ^ 



(^ = l/{ma) is the correlation length in an infinite volume measured in lattice units.) 
The O (a^) lattice artifacts for the step scaling function can be written as 



, , a^ 



u- =u'^ + j^i^iu,z) + 0{a^) , (1.20) 

where the dependence of the correction coefficient z^(u, z) on a/L must be weak (logarith- 
mic). 

It was found in [8] that in terms of this variable the coefficient function u[u,z) for the 
standard action is of the form^ 

usiiu^z) =tQ{u)+ti{u)z. (1.21) 

This is consistent with the In a behavior discussed above. For the non-perturbative con- 
strained (con) and mixed (mix) lattice actions described in section 3 the corresponding 
coefficient functions are given by 

Vconiu, z) = toiu) + tiiu) z + h^u) z^ , (1-22) 

and 

i^miAu,z) = To{u) + Ti(u) z + T2{u) z^ , (1.23) 

respectively. 

The (In of' behavior found for the non-perturbative actions is in apparent contradiction 
with the result (1.12) found for quadratic (perturbative) actions. But one can notice (see 
section 3) the structure 



1 

i(r,.i \ K"" /I °° , 

(1.24) 



8/o(^ioo) V 4 



t2(n) = -8ti(n), 
T2{u) oc ti{u) . 

This structure hints at the fact that here ti{u), t2{u) and 72 (n) are all proportional to 
the matrix element of the same operator and only the behavior of the coefficient function 
(proportional to z for ST and to z^ for con and mix) is different. 

The purpose of this paper is to show by explicitly constructing the large N effective 
action and calculating the relevant matrix elements that the structure described above does 
indeed hold and more generally that Symanzik's effective action description is valid also 
beyond perturbation theory. 



^In [8] it was defined as z = fo{u) + ln(L/a)/(27r). Using eqs. (2.75) and (2.76) one can show that the 
two definitions coincide. Here we want to stress that z depends only on the lattice spacing and not on the 
physical volume. 

*The functions ti{u) together with the similarly defined functions ii{u) and Ti{u) will be given in sec- 
tion 3. 



1.3 Symanzik's strategy 

To verify Symanzik's effective action approach for large N 0{N) sigma models using the 
step scaling function we have to go through the following steps: 

1. Calculate the coefficients in the master formula (1.2) by calculating the 2-point and 
4-point correlation functions in infinite volume on both sides and comparing them. 
In the continuum we need the ordinary correlation functions and also correlation 
functions with dimension 4 operators appearing in the effective action inserted. On 
the lattice side we need the correlation functions expanded near the continuum limit 
up to O (a^) precision. 

2. Calculate the 2-point function (with and without operator insertion) in the continuum 
in finite volume. From this one can calculate the corrections to the mass gap and the 
step scaling function. 

3. Compare the step scaling function directly calculated on the lattice with the ones 
obtained from Symanzik's effective action in step 2 using the coefficients calculated 
in step 1 and verify matching. 

The dimensional regularization calculations necessary for steps 1 and 2 will be pre- 
sented in section 2. The lattice calculations and the matching of the results with Symanzik's 
effective theory will be presented in section 3 and section 4, respectively. 

2 Symanzik theory for the large A^ sigma model 

In this section we give the results for those continuum correlation functions of the 0{N) 
sigma model which are necessary to construct Symanzik's effective action in the large N 
limit. We will use dimensional regularization in 2 — e dimensions. 

We start by recalling the Feynman rules for the large A^ expansion for sigma models 
in Fourier space: 



S"'S (sigma field) propagator: 



■b 



5" 



p^ + m? (2.1) 



AA (auxiliary field) propagator: 



-2 



vAAAAAAAAA/^ 



NB{p) 



(2.2) 



S^S^X vertex: 




a 



b 



6" 



(2.3) 



Here the function B{p) is given by the simple 1-loop integral 

^^^^ " J (2^ (A;2 + m2)[m2 + (A; + p)2]- ^^'^^ 

Finally, correlation functions with r external sigma lines constructed by the above rules 
have to be multiplied by the overall factor 

' f \ '•/2 



and there is an extra 1/2 factor for every closed sigma loop.^ Only sigma loops with at 
least three SS\ vertices need to be considered. 

Next we give the complete list of local operators necessary to construct Symanzik's 
effective action for lattices with cubic symmetry. In addition to the Lorentz scalar operators 
Oi, O2, O3 already defined in subsection 1.1, we need the following ones. 

Us = -^ I -as, (2.6) 

4/0 

U, = ^{I-Sf-^-^, (2.7) 

0/0 0/0 

N 
U5 = —I-I, (2.8) 

Ue = 'i{-j^Os + J2s.df^s}, (2.9) 




M=i 



D 



U7=^-{ -7T^ (Oi + 2O2) + Y, (d^S ■ d^Sf 



(2.10) 



U8 = -^I-S, (2.11) 



N 
Ug = —{d^S ■ d^S) - 2Us. (2.12) 

We will use as our operator basis the operators U3,...,Uq just defined above and the 

combinations 

A^ 1 N 

Ui = -Oi + -[/5 -Us- C/4, f/2 = ^ (Oi - DO2) . (2.13) 

Jo 2 /o 

C/3, C/4 and C/5 are Lorentz-scalar but source-dependent operators, while Uq and Uj are 

components of traceless Lorentz-tensor operators invariant under the discrete (hyper-)cubic 

rotation symmetry group only. As explained in [7], it is necessary to include them in the set 

of operators in the effective action. Here the set is enlarged by the dimension 2 operators 

C/s and f/g. Since their matrix elements are proportional to the dynamically generated mass 

m^, they were not present in the perturbative treatment of [7], but need to be included 

here. 



^Note that loops formed from the legs of our local operators do not count in this rule. All our diagrams 
shown in figures 1-10 are "tree" diagrams from this point of view and only the one corresponding to figure 11 
is a one-loop diagram. The 1/2 rule applies for this last case only. 



2-point functions 

We will use the following notation for connected 2-point functions: 

(5«(x)5"(y)), = g(2)(^,y), (2.14) 

{K,S\x)S'^{y)), = gf\x,y). (2.15) 

Here Ki stands for the integrated (zero momentum) operator 

Ki= I d^zU^iz). (2.16) 

In Fourier space we define^ 

/d^x /"d^ye^^^ie^ygagCs)^^^^) ^ {27,)''5{qi+q2)Gf\k), (2.17) 

('2) 

where gi = — g2 = k and we will also use the "amputated" 2-point functions G) (A:) defined 

by 

An important consequence of the normalization condition S* • 5 = 1 is that 

g^'^\x,x) = l (2.19) 

and this can be rewritten in Fourier space as the gap equation to leading order in 1/A^ 

If d^k 1 

To-Jj2^WT^- ^^-^"^ 

4-point functions 

For connected 4-point functions we will use the analogous notations: 

{S''{x)S-{y)S'{w)S\z)), = g(4) (x, y, w, z) , (2.21) 

{K,S'^{x)S''{y)S'{w)S\z)), = gf (x, y, w, z) , (2.22) 



/d-x/d-./d-./d-.e-e-e-e-a;^)(x,..,.) 



(27r)^ 6{pi +P2 + qi + ^2) Gi ipi,P2, 91, ^2) 



(2.23) 



Mi)( ^ ^ Gr^(Pl,P2,gl,g2) . ,s 

^. lPl,P2,gi,g2j (p2+^2)(p2+^2)(g2+^2)(^2+^2)- ^^'^^^ 



''also similarly without the index i 



Renormalization 

We will have to perforin the usual coupling constant and wave function renormalization, 
which in our problem has to be supplemented by the renormalization of local operators. 

(r) (r) 

We introduce the notation / for the renormalized 't Hooft coupling and G.,L, G)jL for 
renormalized (amputated) r-point functions with and without operator insertion. For the 
coupling renormalization we can use the gap equation (2.20) and we find 

fo = f^'Z{f,e)f, (2.25) 

where 

^ = l + J- + 0{l/N) (2.26) 

and 

f 4-K 27r m 
Here and below the constant 7 is given by 

7 = ln47r + r'(l). (2.28) 

The renormalized coupling depends on the renormalization scale /i, which we will eventually 
identify, for simplicity, with the inverse lattice spacing of the lattice model. 

It turns out that the wave function renormalization constant Z is to leading order 
identical with the one appearing in (2.25) and (2.26) and all our operators Ui renormalize 
diagonally in the large A^ limit. If we define the operator renormalization constants Zi by 

g[2) = Z-'-Z^gM, gJ^Jj) = Z-^/^Z,G^f^ (2.29) 

we have 

Zl = Zq = Z, Z2 = Z4 = Zf = — , ^3 = ^5 = ^6 = ^8 = 1. (2.30) 

Renormalized 2-point functions 

For the construction of Symanik's effective action we will need the dimensionally regular- 
ized, renormalized 2-point functions with and without insertion of the (integrated) local 
operators Ui. The details of the calculation will be presented in appendix A. Here we just 



10 



give a list of the results for the infinite volume renormalized 2-point functions. 

G'g)=/(fc2 + m2), (2.31) 

g[^}^) = 2Tim'^ , (2.32) 

Gg;,)=0, (2.33) 

^3(«) = ^^'(^' + "^')' (2-34) 

G%) = Q, (2.35) 

dX = vik'--Ak')'], (2.37) 



gJ(L)=0> (2.38) 

G^'^. = -47rm2. (2.40) 



Here the notation 



/i=l,2 



is introduced for two-dimensional vectors. 
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Renormalized 4-point functions 

The calculation of the renormalized 4-point functions will be presented in appendix B. Here 
we just list the results of the calculation. Below p = pi + P2 = — (<?! + 92)- 

+ ipiP2){qiq2) - {piqi){p2q2) - (pig2)(P2gi)| , (2.44) 

^^Sr) = -^M^)^P" +Pl + <il + il) , (2-45) 

<l.)-^(^-^ + -^)(^ + -^), (2.46) 

.(4) 



Gg) =0, (2.47) 



^(4) 4/2 f pf ^ p| ^ gf ^ g| 

6(^) A^&(p)\p^ + m2 p^ + m^ g2 _|_ j^2 g2 _|_ j^2 



G 



(4) _8f^( 2( 3 P^ \ I rf/" PiP2 + qiq2 {pip){p2p) + {qip){q2p) 
nm - n\^ \ 4 ^ (^2)2^ +^(, 4 2p^ 



p 
{piP2){qiq2) + {piqi){p2q2) + {piq2){p2qi) 



sr^P^,{Plt^P2^^ + ql^^q2^L)\ , v-^ 

z^ ;;2 j + 2^Pl^^P2^.ql^q2f, 

} , (2.49) 

C) = 1^- (2.51) 

Here we introduced the 2-diniensional limit of the (finite) loop integral (2.4): 

/" d2A: 1 

^P' ~ J (2^ (A;2+m2)[m2 + (A:+p)2]' (2-52) 

and used the notation 

V _ 4Tm^ f 1 at(p) ^ 1 ^ 1 ^ 1 ^ 1 1 _ 

6(p) \ 6(p) 9m2 p2 _|_ ,7j2 p2 _|_ j^2 (j,2 _|_ j^2 ^^2 _|_ ^ J 
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y-_.,„H,,,-4u^'^)m^^i^^i 



4 \ 4 y dm? 32tt \ m? ^ 

2-point function in strip geometry 

After we determine the coefficients in Symanzik's effective action we can calculate the 
lattice correction to the mass gap and the step scaling function from the effective action 
if we know the 2-point functions (with and without operator insertions) in finite volume. 
In this subsection M is denoting the finite volume mass M{L). (The mass parameter 
used in the infinite volume considerations is m = M(oo)). In two dimensions space is 
finite (periodic with period L) in the 1 direction and infinite in the 2 (time) direction 
corresponding to a discrete spectrum of the momentum component ki and continuous 
spectrum for the momentum component k2- In D dimensions we will assume that the 
spectrum of ki is discrete and all other components k2, ■ ■ ■ ,kD are continuous. We checked 
that the D — )■ 2 results remain the same if we assume that in D dimensions only the time 
component /c2 is continuous and all space components k^ (/x = 1,3, ... ,D) are discrete. 
Finite volume quantities will be indicated by an overline. The details of the calculation 
can be found in appendix A. Here we only present the results. 



G% = fik' + M^ 



26(0) 



7(2) _ fh I ,2 ,2 



m 

^2{R) = / ^[h ~ ^2 ~ T 



6(0) 

' + M2) , 
(2) 1 



Gf^n, = ^(k^ + M' 



Here we introduced the finite volume integrals 



2.56) 

2.57) 

2.58) 

2.59) 
2.60) 
2.61) 
2.62) 

2.63) 

2.64) 
2.65) 



h-l ^4^. -d i-f^^^, (2.66) 



(i)(2vr)2g2 + M2 y(L)(2vr)2 {q^ + M^)^ 
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where the meaning of the integration symbol 



(o ^2 (2-67) 

(L) (2vr)^ 

is an integration over the variable §2 and sum over the quantized qi. The above two integrals 
are finite (after regularization) . The parameter p is defined by a regularized (divergent) 
finite volume integral 






(i)(27r)« (g2 + M2)2 D + 2fo 



The finite volume gap equation and the step scaling function 

The gap equation, which follows from the normalization 5" • S" = 1, is of the same form in 
finite volume as (2.20): 

1= M^^ ' - f d^^ ' (269) 

/o J (27r)^ A;2 + m^ 4) (27r)^ k^ + M^' ^ "''^^ 

From this one can calculate [12] the (finite) relation between the infinite volume and finite 
volume mass parameters: 

777/ 

In— + F(n)=0, (2.70) 

where u = LM and 

Fiu)= r^e-^-'(±e-^]. (2.71) 







\k=l 



An alternative representation of the function F{u) is given by [9] 



F{u) = - + lnu-^ + Go(^) , (2.72) 

u \2tt/ 



where 



Go{w) = yi , ^ -I]. (2.73) 

We will often use the function^ foiu) defined by 

in terms of which the step scaling transformation for a general scale s can be written as 

Uu) = Ua{s,u)) + ^^. (2.75) 

For further reference note that /o(tt) has the asymptotic form [8] 

/o(u) = — |-lnu + lnV32| + C'(e~") . (2.76) 



'^The function /o(u) must not be confused with the bare coupHng /o, which has no argument. 
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The finite volume version of the integral (2.52) at zero momentum can also be expressed 
in terms of /q: 



"<»'- 2AP- 
We will need a second function defined similarly to F{u): 



H{u) 



1 



dt 



-tu-' 



2./„ 42 



E 

\k=l 



_k± 

e 4t 



This function has the alternative representation 

■2 1 



U 



H{u) = 7ru + — Inn- -(1 + 2-f)u^ + Att^Gi 



f U \ IT 



V27r 



3 ' 



where 



fc=l '^ 



w^ — k 

2k 



(2.77) 



(2.78) 



(2.79) 



(2.80) 



The rest of the finite volume parameters appearing in the finite volume 2-point func- 
tions (2.56)-(2.65) can be given in terms of these functions: 



uF'{u) 



1 ^^f'(\ 



M2 ^^ ^ M^Hiu) 

h = -—F{u) + ^ 

2-K vru^ 



and 



^ 167r 2 

For the cutoff dependence we will also need the functions /i (n) and /2 (u) [8] : 

,2 



fiiu) 



vr 
6 



12 V27r 



(. 



u 



u u 
12 ~ 16^ 

,2 



, 2^ / u \ 1 



/2(u) = 27rGi ( — )-- + - + —(2k - 1) , 



(2.81) 
(2.82) 

(2.83) 
(2.84) 



where k = -T'{1) -In7r+^ln2. 



Master equation for 2-point function in the strip geometry 

Let us now apply Symanzik's method to calculate the lattice step scaling function by 
writing the master equation (1.2) for the 2-point function in finite volume: 



Qil = y'[s% + a^tu.g% + o{a^)^ 



(2.85) 



From here the lattice correction to the finite volume mass gap is calculated as 

-.2 9 



MLt 



where the coefficients g^ are 






ft = Gl|)(0,iM). 



(2.^ 



(2.87) 
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If we take // = 1/a, i.e. identify the renormalization scale /x with the inverse lattice spacing, 
we can give the connection between the renornialized 't Hooft coupling / and the lattice 
parameters using (2.27): 

2,71 ^ 

lne = y-^, (2.88) 

where ^ is the lattice correlation length in infinite volume. 

Using the results of the previous subsection, the coefficients in (2.86) are 

^^ 26(0) ^ ^ 

52 = 4^|=|t-m4, (2.90) 



2 U(0) 
53 = 54 = 55 = 98 = , (2.91) 

g, = 2fM' - 1^ , (2.92) 

57 = -252, (2.93) 

Oq = -= — . (2.94) 

^^ 6(0) ^ ^ 

3 Lattice actions 

We consider here the large A'^ limit of the standard, the constrained and the mixed lattice 
actions. Since the standard and the constrained actions can be obtained as special cases 
of the mixed action, we shall discuss in detail only the latter one. 

3.1 Standard action 

The action is given by 



At[S] = ^5](5^S.)2, (3.1) 



2 

where S^ = 1 and d^ denotes the forward lattice derivative. 
As usual, the large N limit is taken as 

/3st = ^, /st = fixed. (3.2) 

Jst 

Note for further reference that for large A^ the distribution of {d^Sx)'^ approaches a 6- 
function at (d^Sx)'^ = /st/2 (for D = 2 Euclidean dimensions), as can be shown in pertur- 
bation theory. 

For the standard action the infinite volume mass gap m is given by [12] 



fst 

where z is given by eq. (1.19). 



l + 0{m'^a'^) , (3.3) 
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3.2 Constrained action 

The action is the same as the standard one except that the configurations are restricted 
by the constraint (d^Sx)'^ < e: 

[ oo otherwise , 

with /3con = N/fcon- As has been shown in [8], the constraint (S^S^)^ < e is effectively 
replaced by the (^-function (5((9^Sa;)^ — e) in leading order in 1/A^. As a consequence, the 
physics for e < /con/2 (to this order) is equivalent to the topological action [8] where only 
the constraint is present, i.e. 



.4„,[S,.|° '-(W^<- ,3,5) 



cxD otherwise . 



The mass gap in this case is also given by (3.3) with fst replaced by 2e. (For e > /con/2 the 
constraint is irrelevant in leading order, and the physics is given by the standard action, 
with /st = /con- We assume here that e < /con/2.) 

3.3 Mixed action 

The mixed action [8] is given by 

« ^, 

(3.6) 



AmiAS] = ^)^{d,S, 


^f + ^r i^iid.sxf? 


X,fJ, 


X,fl 


The large A'" limit is taken as 




N 


2N 


Pm — r ) 
/m 


7m — 9 1 



(3.7) 



keeping /m and k^ constant. 

The infinite volume mass gap is determined by the effective coupling /m defined by 



11/11 , , 



/m 2/^ V 4/ 



m "'m 



and is given by the same expression as (3.3) with /st replaced by /m- 
Besides /m we introduce 

rra = -p- and Qm = -^ (rm + V'>^m + 4] ■ (3.9) 

Jm -^ ^ ' 

While the mass gap is determined by /m alone, the cutoff effects depend on the ratio r^ 
as well. 

It is easy to see that rm = oo (gm = oo) corresponds to the standard action. The 
choice Tm = {q^a = 1) gives the "purely quartic" action (where /Sm = 0). It is less obvious 
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that the r^ — )■ — cxd (gm — ^ 0) hmit gives the constrained action. For large negative rm the 
action density in (3.6) has a deep minimum at 

(5^S.)2 = -^ = IfJl + O (r-2) ) , (3.10) 

which shows that this hmit indeed corresponds to the constrained action with e = /m/2. 

The technical details (introduction of auxiliary variables, the gap equation, the calcu- 
lation of the 2-point and 4-point functions) can be found in [8], and we shall not repeat 
them here. 

The cutoff dependence of the step scaling function S(2, n, a/L) is described by (1.20), 
(1.23). The functions Tj(ii) are given by 



Ti{n) = j^^ (^Mu) - \'^^{u'^)^ , (3.11) 



ff. f \ f / \ , ^ "^f f \ 2/2 (n -u fo{u) / w. X , 1 
^o{u) = fi{u) + -u fo{u) 2—— u/o(n) + — 



i + q^ V ""^ ' 47r, 

^^(^u) = --u', (3-12) 

8 

The functions fi{u) are defined in (2.74), (2.83), (2.84). Note that the coefficients Ti{u) 
depend on the parameter q^ defined in (3.9). As shown above, the coefficients ti{u) and 
ti{u), describing the cutoff dependence for the standard and constrained actions can be 
obtained from Ti{u) by setting q^ = 00 and q-m = 0, respectively. For the standard action 
(gm = 00) the z^ term is absent. 

4 Matching to the lattice results 

As described in subsection 1.3, we calculate the 2-point and 4-point correlation functions 
on the lattice in infinite volume up to O (a^) , and compare these with the results obtained 
using Symanzik's effective action (1.1). 

4.1 Matching the infinite volume 2-point function 

In leading order in 1/A^ the 2-point function for the mixed action is the free lattice propa- 
gator 

where ttt-q is expressed by the infinite-volume mass gap as 



2 , 'ma- , . N 

mo = - smh — - , (4.2) 



and F is related to the effective coupling /m by 



L = F r^^ + ^rl + A_Aa^^ ^ ^ ^ ^ ^^,^,^ ^ ^ ^ ^ (exp(-4^/F)) . (4.3) 
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The gap equation is 

where 2 is defined in (1.19) and 

We now have 

G(fi(A;) = -(A;2 + m^) + — |t^(A;^ - m^) - ^m\k'^ + m'^)\ + O (a^) . (4.6) 

According to Symanzik's conjecture this should be equal to 



/ 



I Gg) (A;) +a'Yl ^«^1|) (^) + ^ («') | • (4-7) 



Using (2.31) we find that the scaling parts match if the wave function renormalization 
constant is of the form 

y2 = i- (1 + wsa^ml) + O (a'^mt) . (4.8) 

Matching the O (o'^) terms gives 

^(A;^ - m^) - i;n?{k^ + n?) = wsm^k^ + m^) + i ^^.gJ^]j)(A;) . (4.9) 

Using the formulas (2.32)-(2.40) we find 

(4.10) 
(4.11) 
(4.12) 
(4.13) 

4.2 Matching the infinite volume 4-point function (standard action) 

We analyse first the ST case. In this case the 4-point function is 



V6 = 


1 
24' 




V3 , 


1 
16' 




, ^8 
"^ + 2m2 = 


16 ^ 


-m, 


2V9 

Vi 2 = 


/ 
967r' 





NBi^ttip) pI + mlpl + ml qf + mg g| + mg ' 



C^*^ = fJi^t - t t t (AU) 

='latt A7D (-^A .;2 , ,„,2 .;2 , „.2 ~2 , „„2 ,-2 , „„2 ' V^--^^^ 



where p = Pi + P2 = —Qi — Q2 and 



7r/a j2 1 

Aatt(p) = / 7^ r ^^^ • (4-15) 



n/a (271") ^^2 _^ ^2-) ^^2 _|_ (^TlfT^) 
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Expanding in the lattice spacing we have 

BiMp) = Kp) + aHiip) + O {a'') . 

The functions b{p) and bi{p) are given in appendix C. 

The amputated 4-point function on the lattice is then 



(4.16) 



G 



(4) 
latt 



+ 



2a^ 



Nz%{p) Nz'^b{p) 



''-'-"-§ 



1 
12 



4 44 44 44 4n 

PI — m P2 — Tn ql — m q^ — rn 
p\ + m? P2 + m'^ qf + m'^ g^ + "^^ 



(4.17) 



+ 0(a' 



which we have to match to the DR result 
1 



22 v^ + '^wso^ml + O (a^TTT-Q 



Z2/ 



))|Gg) + a^5:.4^))+0(a^)| 



(4.18) 



Using (2.42) we see that the scaling pieces indeed agree and from the rest we have 
1 



2^-' + ¥4 

b{p) 



12 



4 44 44 44 4n 

p| — m P2 — m ql — m g^ ~ "*- 



Nb{p) 



- 2w;8m" + ^^^ Yl ^*^S) • 



(4.19) 



Here we have to use the results (2.43)-(2.51). Because contributions of some operators have 
momentum dependence not occurring in the above expression we first have 



Next we get 



V2 = V4^ = Vt = 0. 



1 



^3 = 4 



which combined with the results of the previous subsection gives 

f5 = and Da = —2rr' 

16 

The last matching equation is 



and ^8 = — 2m^(^ + wg) ■ 



vi m^dbip) 1 



(4.20) 
(4.21) 

(4.22) 
(4.23) 



After inserting the expressions for b{p) and bi{p) given in appendix C, we obtain the simple 

f(. i_Ul + ifl„8_i_r'(l)V (424) 

Eqs. (4.10), (4.20), (4.24) are consistent with (3.30) from [7] where we found to first 
order in coupling constant PT^: 



"^=4^ 



V^ = -^ + OU), ^6 = ^ + 0(/), VT = 0{f). 



(4.25) 



Note the extra A'^ factor in our U\ with respect to f/i of that paper. 
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4.3 Matching the infinite volume 4-point function (mixed action) 

Referring to appendix B of [8] the mixed case amounts to adding an extra piece to the 
propagator of the auxiliary field, 

A{p) ^ A(p) + a2 J^X^(pi,p2)X.(gi,'Z2)A^.(p) 



H,!/ 



where 



A(p) 



^^.AP) 



1 



Blattip) ' 

2tt " 



F2 



h 6, 



1 



1 

k2 = 



+ k2 

1 



+ 0{a^) , 



and 



TT - 2 + TTq^ ' 

X^,{pl,P2) = -pipiP2pi 



27r(l + ql, 

hf,{p) 
h{p) 



(4.26) 

(4.27) 
(4.28) 
(4.29) 

(4.30) 



Here h{p) = -Biatt(p) and hfj_{p) are given in (4.15), (C.31). 

Expressing the terms with specific momentum dependence by the appropriate combi- 
nations of G ,L. , after some algebra we get for the additional part of the 4-point function 
of the mixed action 



Y,X,{pi,P2)X,{qi,q2)A,,{p) = ^{^ (GJt«) - 2^: 



fl,!/ 



.(4) 



^'='7 






^ "(4) ^ 1 ^(4) ^2^2 (4) 



+ ^G 7.^ + ^G 7.^ + 



-Gl 



The matching gives 



f2^3(R) ^ p^MR) ^ f2 ^8(R) 
;f^ = vf+Xi, i/5. 



)} 



+ 0(a' 



where the Xj are determined through 

Z^ ^i^i{R) - ^ 



(4.31) 



(4.32) 



Y,X^{p^,P2)X,{qi,q2)A^u{p) + 0{a^) , (4.33) 
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and so we get 



<" = / 



-A7rz^k2 H 

4 327r 



(4.34) 



^^2-^'^ = ^ , (4.35) 



I 

4 



^mix ^1 4^^^^ ^ (4_3g) 



<^ = -^, (4.37) 

<" = -4 + 2^^^2, (4.38) 
lb 

vr = ^ , (4.39) 

vf- = -^ , (4.40) 

vf"" = -2m2 [V; + 2iTzk2 + ws] , (4.41) 



2 
777 
mix 

«9 -^ 






(4.42) 



The coefficient v^^^ can be determined only up to the parameter ws which can be fixed if 
we first impose a renormahzation condition on the lattice 2-point functions. 

Note that the matching coefficients for the mixed action are non-perturbative in the 
coupling constant, some of them contain a 1// factor. This explains why the leading artifact 
in this case is a^ In a as opposed to a^ In a for the standard action. The mixed action is 
perturbative (i.e. the quadratic part dominates) only when /m ^ Km, i.e. Qm S> 1. In 
this case the q^ in the denominator of (4.29) compensates the 1// factor in the continuum 
limit. 

To minimize the lattice artifacts for the step scaling function one has to take [8] q^ = 
8//m + (1). Inserting this and z = 1// + (1), and l//m = l/Z + O (1) into (4.34)-(4.42) 
the coefficients vi and v^ become O (/) while for the standard action they are O (1). As 
has been shown in [7], indeed the operator Ui is responsible for the leading lattice artifacts. 

4.4 The matching for the step scaling function 

Let us parametrize the "couplings" gi defined in (2.89)-(2.94) as 

-g, = 2rM\L)4jTr-y (4.43) 
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Using also the finite volume results from appendix B of [8], we have 

ei = 0, ^>^{u) = -—, (4.44) 

62 = 2, ^>2{u) = -1 Ti- + uf^{u)\ {2h{u) - u^Uu)) , (4.45) 

66 = 1, ^e{u) = 24/i(n) + 3u^fo{u) - -— , (4.46) 

ovr 

67 = 2, ^j{u) = -2^2{u) , (4.47) 
69 = 0, *9(w) = L^ , (4.48) 



and 



*3(i() = ^i{u) = ^5{u) = ^siu) = 0. (4.49) 



It is interesting to observe that the off-shell operators U3, f/4, C/5 and Us do not contribute 
to the finite volume mass. Ug does contribute, but since 'ifQ{u'^) = (2L)^, its contribution 
to the step scaling function also vanishes: 

*9(n) - \Mu'oo) = 0. (4.50) 

Using the mass formula (2.86) we find 

with 

■^(u) = y^i(n) + f{v2 - 2v7)^2{u) + ve^eiu). (4.52) 

Inserting the coefficients Vi given by eqs. (4.34)-(4.42) one obtains for the mixed action^ 

2-Tr 1 

V2-2vr = ^—^7, (4.54) 

vr - 2 + 7r< / 

^^6 = ^ . (4.55) 

These results agree with those of [8]. Note that matching the finite volume mass gap does 
not fix completely the coefficients Vi ~ it yields only the relations given above. Also note 
that the coefficients of the effective action obtained by matching in finite volume do not 
depend on the volume, as expected by general considerations [13, 14]. 



®The results for the standard and constrained actions can be obtained from here be setting qf^ to 00 and 
0, respectively. 
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5 Conclusions 

The extrapolation of lattice data to the continuum limit is a crucial step in the deter- 
mination of low energy physical quantities in Quantum Chromodynamics. Our present 
theoretical understanding of lattice artifacts is based on considerations of Symanzik. He 
argued that asymptotically the leading ultraviolet cutoff effects are described by an ef- 
fective continuum Lagrangian containing a finite number of higher dimensional operators, 
restricted by symmetries only of the underlying lattice regularization, and with coefficients 
depending on the lattice spacing. 

In this paper we have given additional support for Symanzik's theory of lattice arti- 
facts, by determining the local effective Lagrangian for a class of lattice actions of the 
2-dimensional nonlinear O(A^) sigma model in a non-perturbative setting of the 1/A^- 
expansion at the leading order. The class of models considered includes also ones which 
are considered to belong to the same universality class as standard actions, but for which 
the usual perturbative expansions are (as yet) not available. 

The effective actions are shown to reproduce previously computed lattice artifacts of 
the associated step scaling functions which are defined in finite volume. Once established 
the effective action can be used to predict cutoff effects for various observables. 
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A Two-point functions with operator insertion in finite volume using 
dimensional regularization 

Ordinary 2-point function and source-dependent operators 

From the leading order large N Feynman rules we see that the sigma field propagator is of 
the free form 

G^^^ = fo{k' + M^). (A.l) 

We see that the wave function renormalization constant Z is the same as appearing in the 
coupling renormalization and the renormalized 2-point function is 

G% = fik^ + M^). (A.2) 

The correlators of the source-dependent operators C/3, C/4 and C/5 are given in [7]. The 
correlators of U2 are simply related to the ordinary correlation functions (without operator 
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insertion): 

Gi''\p^, . . . ,pr) = \(J2 pI) ^~^'^(P1' • • • 'Pr)- (A-3) 

From this representation it is clear that we have for the operator renormalization constant 
Zs = l. 

The correlators of C/4 are related [7] to those of the isospin tensor operator S°'S'^ — 
-^5"^. Its 2-point correlation function vanishes and therefore the operator renormalization 
constant Z4 can only be calculated from its 4-point function. 

The operator C/5 corresponds to a contact term: its 2-point function is given by 

G'i^=Mk^ + M'f (A.4) 

— (r) 

and all higher correlation functions vanish: Gg =0 for r > 2. We see that Z5 = 1 and 



le renormalized 2-point function is 






Gf^R^ = fik' + M^f. 




(A.5) 


inally for the source-dependent operator C/g we have 






g(^) = Jg» 




(A.6) 


)r all r. Prom the r = 2 case we obtain 






G^' = ^-^{k- + M% Zs = l, Gf,^y- 


= {{k' + M^). 


(A.7) 



The coupling-related operator C/g 

This operator is related to the integral of the action density 

-L fd^xd^s-d^s (A.8) 

and if we denote the corresponding correlation functions by Q^ , where X is any combination 
of operators, we can derive the following Ward identity: 



^^'"''W.^'-fW^^"- <"' 



Applying this to the 2-point function we get 



^(2) fo_ _ 

k'^ + M^ B(0)(fc2+M2)2' 



O'a' = tt'^. - ^7:w7T:^r-^j^^ (A-io) 



where we used the relation 

9M2 1 

h-j^ = T^TTT, (A.ll) 

which can be obtained from the gap equation (2.69). Here B(0) denotes the finite volume 
version of the integral (2.4) evaluated at p = 0: 
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Ol; 




Figure 1. 

A contribution of the 2-point function with insertion of the operator Oi or ©2- The two 
black blobs denote the two signia- field contractions in (1.5), and have similar meaning in 
all other figures. The isospin indices of the outgoing legs are contracted. 



Ol: 



^ 



Figure 2. 



Using the definition (2.12) we find 



From this we see that Zg = Z and 



G\ 



(2) _ 1 



5(0) 



^9(«)- 5(0)- 



(A.13) 



(A.14) 



The operators ?7i, U2 

Let us introduce the index notation 1, 2 for correlation functions with insertion of the space 
integral of the local operators 

1 






i = 1,2. 
We then have for the 2-point functions (see figures 1, 2): 



(A.15) 



(2) 



G' 



77(2) 



f2 / ^ 






5(0) 



i/=i 



C, 



B{0) 



(A.16) 
(A.17) 
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where 

r d^q ql r d^q q^ 

These integrals can be expressed in terms of the differences 

ni-'H2 = h, Ci-C2 = ( (A.19) 

and sums 

^ M2 ^ 1 _ 

E^M = -^' J;/:^ = ^-M2s(0). (A.20) 

Taking into account that there is just one momentum component (gi) with discrete spec- 
trum and D — 1 continuous ones, we have 

D 

Y.'^t^ = Dn2 + h (A.21) 

and similarly for £^. 

Using the above relations we can now write 

Gf = -^^(k^-^j:C.]=-^e + M^) + ^ (A.22) 



leading to (using (2.13)) 



^(2) _ ^' 7 _ 7 77(2) _ M2 

2my ^^-^' ^^(^)-^)- 



&{^ = -^^^^, Z^ = Z, G\'>^-^ = ^^^^. (A.23) 



— (2) 

An alternative way to express G^. is 



1 2 1 ' ^ 5(0) B(0) B(0)J ^ ^ 

Comparing to 



2 1 ' ^ 5(0) 5(0) 5(0) 



Gf = 4^ { n^k' + /^fc? - ^^^^ - '"^fJ::'' - 4^ \ (A.25) 



we obtain using the definition in (2.13) 

Gf = ^{k'- Dkl + {D- 1)=^} . (A.26) 

Here we can demonstrate that our results are independent of the way we regularize the 
finite volume problem. Indeed, if, instead of (A.21) we write 

D 

E ^^ = Dn2 + {D- l)h (A.27) 

^=1 
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Figure 3. 



corresponding to D— 1 discrete and one continuous niomentuni components, we find instead 
of (A.26) 

G?) = ^ S^Dkj -k' + iD- 1)=^} . (A.28) 



In the D -> 2 limit, both (A.26) and (A.28) lead to 



Z2 = l, G?^m = ^\kl-kl + =^\. (A.29) 



The integrals (2.66) corresponding to the differences h and i are finite. 

The tensor operators Uq, U-j 

The main part of U^ is the term 

D 



^ V {d^S ■ d^Sf (A.30) 

■^0 ^ 



and corresponds to the 2-point function 

4/o^E^.(^^^)=8Gf. (A.31) 

Combining this with the remaining contributions from (2.10) we find 

G? = f 8 - -^] Gf - -^ Gf = -^G?^ (A.32) 



and hence 



1 7^(2) _ o7^(2) _ ^2, J ,2 ,2 



^7 = ^, G^(k) = -2G-^) = f^hJ^kt-ki-^^}. (A.33) 



The main part of Uq is 



/V ^ 



and corresponds to the 2-point function (see figures 3,4) 



^H' 



2/o|i:^;^-ij^E^.|' (A-35) 
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Qi 



• \aaaa/^<^ 



Figure 4. 



where 

Using the definition (2.9) the 2-point function of Uq becomes 

Gf = Vo^K + ^ {M' - (kr) -^% (A.37) 

with 

^ = E^M + 7y-:^^=P + 0(e). (A.38) 

t^i D + 2 fo 

The final result for this operator is 

Ze = 1, 4(k = Vk' + ^(M^ - (fc^)^) - Ig. (A.39) 

To obtain the infinite volume results (2.31)-(2.40) from the above formulae is straight- 
forward. One has to use the L — )• oo limits 

1 3m? 
6(0) = - ^ and p= . (A.40) 

Calculation of the finite volume integrals 

By introducing a Feynman parameter, we can rewrite the gap equation 






Jo Jo J{L) (^^j 



and similarly the basic loop integral at zero momentum: 

/■oo I- aD , 

B{0)= I tdte-'^^' I T^e-'^'. (A.42) 



Let us concentrate on the contribution of a single one dimensional integral in the last 
D-dimensional integral: 

ri!5e-'.= = '. (A.43) 

J_^ 27r Vi^ 

The corresponding "integral", if the momentum variable is discrete, is really an infinite 
sum of the form 

^ oo 
n=— oo 
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Here and in all other similar sums we encounter the theta-function S{x) defined by 

oo 

S{x)= Y^ e""^"'. (A.45) 



Using this notation, (A. 44) can be written as 

1 fAirt 

A useful identity satisfied by the function S{x) is 



(A.46) 



(A.47) 



which can be proven using the Poisson resummation formula. Thus an alternative form of 
(A.44) is 

^ S(^]. (A.48) 



'^^^-^x'il 



Now we rewrite the gap equation and the zero momentum loop integral as 

/o Jo 
BiO) 



(^\ 



(4^t)-D/2 \4TTtJ 
f2 1 



dte-^^'—^TTT^. 
(47rt)^/2 ' 



tdte"*^ 

(4^t)^/2 ^ V4^iy' ■ 



s(.'' 



(A.49) 
(A.50) 



The advantage of using the infinite sum S{x) with this argument is that the n = term 
is identical with the infinite volume limit of the given quantity, while the n ^ terms are 
UV finite [14]. 

From (A.49) we obtain the equation 



dt e 



-tM'^ 



-tm^ 



1 



+ 



(47rt)^/2 7o 



die 



^tM'2 



1 



(47rt)^/2 



A-RtJ 



0. (A.51) 



Both integrals have finite Z) — )• 2 limits and evaluating the first one and using the definition 
(2.71) in the second one leads to the final result (2.70) for calculating the LWW coupling 
u = LM. The expression (A.50) is also finite for D — )• 2 and gives 



AnM^ 6(0) =u'^ die-*"" < 1 + 2 ^ e"^ f = ^ ~ "^'M- 



Next we consider integrals of the form 
d^(? ql 1 



fc=i 



rdte 



-tM^ 



d^ 2 -tq^ 



i^L) (2vr)« ((?2 + M2)-+i r{a + 1) J, ^ ^^ J^^) (2vr) 

We will need integrals with cr = 0,l. If/^T^l then the last integral is simply 



1 



1 



Si^] 



2t (47rt)^/2 \i7rtj 



(A.52) 



(A.53) 



(A.54) 
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whereas for // = 1 it is 



1 



S 



L' 



's(i^). 



2t (4^t)^/2 V4vrt/ (4^t)-f5/2 dt \4TTtJ 



(A.55) 



To evaluate h and I defined in (2.66) we only need to integrate the difference between the 
above two terms, which is finite as -D — ?■ 2: 



1 d 



Airt dt 

2^ Jo 



V47rt 



1 



die-*"' {t'^-^n^ + (1 - a)r-^} ^ e" 



(A.56) 



k=l 



For cr = 1 we obtain from this the first and for a = the second relation in (2.81) 
Finally we consider 



tdte 



(L) 



(2vr) 



9.6 



-tg^ 



D '^l^ 



Here for /-i 7^ 1 the last integral is 



S{ — 



4*2 (47rt)^/2 \A7Tt 



whereas for /x = 1 it is 
1 



3,/L^ 



1 a /L^\ a^ /L^\ 
(47rt)^/2 [4t2 ^ y4TTtJ t dt \4TrtJ dt'^ \4TrtJ 

Summing over fi, as required by the first term in (2.68), we get 



tdte 



'tM- 



1 



SD^/L^ 



1 a /L^\ 5^ /L^ 
(47rt)^/2 [4t2 ^ \47rtJ t dt \4TTtJ dt^ yint 



This has to be combined with the second term in (2.68), which is of the form 

D + 2J0 ^^ (47rt)^/2 "^ V4vrty' • 



(A.57) 



(A.58) 



(A.59) 



(A.60) 



(A.61) 



Let us now separate the contributions into a "finite" part, which is obtained by the re- 
placement 

(A.62) 



'[^t) 



^'isj-' 



corresponding to UV finite integrals and an "infinite" part, which corresponds to the sub- 
stitution 

■ 1. (A.63) 



'[17,) 



We first evaluate the "infinite" contributions. We find 
3^ f^ , tAf2 1 n/2 6M2 '■°° 



4(47 



(Z) + 2)(47r)«/2 7o 



dte-**^'t-^/2 



^{fr,-«/..^ra-«/. 



3M 



D 



(A.64) 



(L> + 2)(47r)^/2 



T{2-D/2). 
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In the D — 7- 2 limit this is finite and gives SM^/lGvr. The "finite" part gives 



tdte 



^tAP 1 



A-Kt 1 2*2 



ATTtJ 



1 a /L^\ 9^ /L^ 

tdt \4TTtJ df^ \4TTt 



+ 



3M^ 



die" 



-tM^ 







M2 

M2 



die 



-tM^ 



1 

ivrt 

/ CO 

E 



*'S)-' 



1 



^"^^^ 1^+^+"' 






-M^ 



h 
2' 



Adding the two contributions we finally have the result (2.82). 



(A.65) 



B Four-point functions with operator insertion in infinite volume using 
dimensional regularization 

Ordinary 4-point function and source-dependent operators 

The leading large N 4-point functions are of order 1/A^. Using the Feynman rules given in 
section 2, the leading order ordinary 4-point function is 



^(4) 



2/g 
NB{p)' 



(B.l) 



After renormalization we get for the renormalized 4-point function (2.42). The 4-point 
function of the source-dependent operator C/3 can be obtained from the above by multiply- 
ing it by a momentum-dependent factor: 



G 



(4) 



^° (pI+pI + qI + qI) 



2NB{p) 



(B.2) 



and after renormalization this becomes (2.45). Correlation functions of the source-dependent 
operator f/4 can be calculated using the rules given in [7]. For the 4-point function we get 



Gf = !l(^ + rrA('i±^^rn- 



(B.3) 



This allows us to determine Z4 (the only operator renormalization constant not determined 
by the 2-point functions). We obtain Z4 = 1/Z and (2.46). Next we recall that 



^(4) _ ^(4) _ r. 



(B.4) 



since all r > 2 correlation functions of U^ vanish. Finally the 4-point correlators of C/g 
coincide with the ordinary ones: 



G 



(4) 



NB{py 



G 



(4) 
8(R) 



2P 
Nb{p) ' 



(B.5) 
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Figure 5. 



The action-related operator Ug 

From the Ward identity (A. 9) we obtain 

where the coefficient X^°' is the bare version of (2.54): 

X(°) = ^ r^95(p) 1 1 , 1 , 1 ] ^ ^ , O .,^ 

Bip)B{0) \ B{p) dm^ ^ pI + 7v? pI+7v? ql + m'^ ql+w?\ ^ ' ' 

(B.7) 
Subtracting the U% part we get 

and after renormahzation (2.51). 

The operators Ui, U2 

There are two types of contributions to the 4-point functions G- and G- (recah these 
are correlators of the operators (A. 15)): either the four external legs are attached to the 
four legs of the local operators, or all four external legs are coupled to one pair of operator 
legs and the other pair of operator legs are contracted with each other. The latter type of 
contributions are proportional to the action-related 4-point function G^ . Indicating the 
first type of contributions by a hat, we have 



(^(4) ^(5(4) _ ^^(4)^ ^(4) ^(^(4) _ rn^^(4)^ ^^, 



The "hatted" contributions are constructed from the building blocks shown in figures 5-7 



^i'^ = m ^A'mCpI'^s) T,,{q,,q2), (B.IO) 

'^? = ^T^-iPl^p2)T^,{qi,q2), (B.ll) 

where 

2 

-Tfj,y{pi,p2) = PitxP2v + PivP2fi + -j^rriM^tJiv + BoPi,p„). (B.12) 
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1,2 



P 

AAAA/^< 




-Pi 

•P2 

■^1 
92 



Figure 6. 



o 



1,2 



P 



\JVv\AA/^- 





Pi 

P2 
l^l 
92 



Figure 7. 



The p-dependent scalars ^Qj ^o are defined by the regularized loop integral 

[ d^g qtJ,{Q+p)u _ , , „ 

J {2tt)D {q2+rn^)[^^ + {q+py] "^oV + ^oP^P. 

and are explicitly given as 

1 m'^+p'^/4 



An 



P^Bo 



2{D - l)/o D - 1 

D-2 



B{p), 



+ 



2{D - l)/o 
^0 is divergent but Bq has a finite Z) — )■ 2 limit: 

Bo = fco + O (£) , 
Using 



m^ p^{2-D) 



D-1 4(L>- 1) 



B(p). 



p^6o = -1 1- m^Hp). 

47r 



1 



D^o + p' Bo = ^ - (m^ + pV2) B{p) , 

JO 



=.(4) 



G- can be simplified: 



r(4) -ilL,2, Pi+P2 1 



2 , gl+gj 1 



/oi?(p) 



/o 2m^/o^ ^ ^(4) ^ ^(4) 



NB^{p) NB{p) 



(B.13) 



(B.14) 



(B.15) 
(B.16) 



(B.17) 
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Comparing the last line and the definition in (2.13) we see that the 4-point function of Ui 
simplifies to 

(4) _ /o "i^/o 



G'i 



-X^°\ 



(B.18) 



NB'^{p) N 

The renormalized version is obtained by the substitution fo^f and replacing the rest by 
its finite Z? — t- 2 limit as given in (2.43). 

There is an alternative representation of G- which will be useful later. Going back 
to the definition (B.IO) we can write 



G 



(4) 



/3 (AD'^Al , 4DAo 



+ 



4iV [ B2(p) B{p) 



piP2 + qiq2 + 



2Bop^ 
B{p) 



+Ti 



where 



-Fi 



iBiip'f ^ ABop^ 



1 B^ip) ' B{p) 

Similarly, from (B.ll) we calculate 



{piP2 + qiq2)+ HpiP2){qiq2)- 



^(4) /g fiPAl 4Ao 
2 An\b^{p) B{p) 



piP2 + qiq2 + 



'2Bop^ 
B{p) 



+ F^ 



(B.19) 
(B.20) 

(B.21) 



with 



-^2 = ^^BHp) + 1^ [(PiP)(:P2P) + {qip){q2p)] + 2{piqi){p2q2) + 2{piq2){p2qi). (B.22) 

To compute the 4-point function of U2 we first observe that the action-related terms cancel 
and we can start from 



G 



(4) 



G 



(4) 



DG^ 



(4) 



(B.23) 



'2 "1 -"2 

After some algebra we see that the divergent j4o-dependent terms also cancel and we are 
left with 



G 



(4) _ /o 



N 



D 



(„-i)i«?4V+^°''' 



p. 



B{p)J B{p)V 
D 



D 

P1P2 + qiq2 — n 

pZ 

D 



{PIP){P2P) 



(B.24) 



-^{qip){q2P) + {piP2){qiq2) - TriPiqi)iP2q2) - -^{piq2){p2qi 



The renormalized 4-point function (2.44) is obtained from this in the L> ^- 2 limit by 
making the substitutions /o — >• / and 



Bop" 



hop^ 



B{p) b{p) 

The tensor operator C/7 

The essential part of Uj is the operator 

iV ^ 



m 



47rb{p) 



Q. 



(B.25) 



- Y,id,S ■ d,Sf 



h 



(B.26) 



M=i 
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and we will denote its 4-point function by G^ . Very similarly to the previously discussed 
cases, we can separate this 4-point function into two parts: 

Gf=Gf-'-^G^i\ (B.27) 



where 



^7 - ^ Z. [P^^P^^ +B^) + B(p)) ['^'^'^'^ ^W) Bip)) • ^ ^ ^ 



This can be simplified to 

2Sop2 



^(4) 8/3 f DAI I ^0 
7 AT \b^{p) B{p) 

where 



PiV2 + q\q2 + 



B{p) 



+ F^ \ (B.29) 



Using (2.10) we notice that action-dependent terms again cancel and we can write 

^7 -G^ -;DT2^i ":dT2^2 • (B.31) 

The divergent, ^o-dependent terms also cancel here and we are left with 

(4) ^ fl f _ ^_jr _ -^T.) . (B.32) 

The final result (2.49) is obtained from this last formula by making the replacement /o — ?■ / 
and writing the D — t- 2 limit of the rest as a quadratic polynomial in Q. 

The tensor operator Uq 

Here we will use the following representation [7]^*^ of the operator Uq: 

6 / 47V \ 

U6 = Qe- -^^ (U^ - 4C/4 + —OA . (B.33) 

Here Qq, defined in (A. 34), is the main part whose correlation functions we will denote 
by the dotted index 6. Using the fact that U^ is a pure contact term from the above 
representation we have 

0«*' = 4*'-I^(or-<)- (B-34) 

There are several contributions to the four-point correlation function G- . The corre- 

D 

spending Feynman diagrams are shown in figures 8-11. 



"This follows from formulae (A. 10) and (A. 34) of that paper in the large A^ limit. 
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Figure 8. 




Figure 9. 




Figure 10. 




Figure 11. 



The sum of the contributions of the Feynman diagrams shown in figure 8 is 



4/0 V^ 



4 4 



4 4 



NB{p) ^ I p^ + 771^ P2 + ^-^ 9l + '^ Q2 + "^ 

/i — 1 v. 



(B.35) 
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The 9 type contribution is 



^f^f^^y 4 (B36) 

NB-^ip) J (27r)^ ^ (g2 + m2)2 [m^ + (g + p)2] ' ^ -^"^ 

/J, — i 

The sum of the 10 type contributions is 

iVS(p)S(0) \^f +m2 pl+m"^ ql + 'm? q2+m'^i_.J {2'k)'^ {q^ + m?Y 

(B.37) 
and finally the 11 type contribution is 



NB'^{p)B{0) dm? \^ J {2tt)D {q'^ + m 
\f^ — ^ 

Putting together all contributions we have 

^ rO D ( 4 4 4 4 

^(4) ^ 4/2 ^ ( Pi^ ^2^ glM , g2^ 

^6 A^Rfr)^ Z^ 1 Ti2_Lm2 "^ „2,^2 "^ „2_l™2 "^ 



NB{p) ^ I p2 _|_ ,7^2 p2 _|_ ^2 g2 _|_ j^2 (^^m? \ 

+ MxW T [ ^^ ^J fB 391 



/VR2/'„A Z^ 



iV52(p) Z^ y (2vr)S (g2 + ^2)2 [^2 + (^ + p)2] " 

Using our earlier results we have 






and finally 



A ^0 D ( 4 4 4 4 

^.(4) ^ 4/g y^ ( Pi;. ^ ^2^ ^ glM ^ g2M 

+ (D+'2TL(p) (^^+^^ + ^^ + ^2') 

,4/g^^ro^ff / d^, .^ 6 m^l (^-^1) 

^ 1^7 (27r)^ (g2 + m2)2^Z? + 2 /o 

, 8/0^ [^ f d^, g;^ 3^ 

iVS2(p) ] 4^ J {2n)D (g2 + ^2)2 [m2 + (g + p)2] (D + 2)/o 

The coefficient of X^°' can be simplified if we evaluate 

J (2^ (q^+m^y " *(^a/3'^7'5 + '^"7'^/3<5 + ^aS^fj^). (B.42) 
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^ can be calculated by considering the contraction 



This gives 






(2vr)^ A. (g2 + ^2)2 (D + 2)/o ^ I) + 2 ■ ^ • ^ 

/J, — 1 

The four-point function is now simplified to 

® NB{p) ^ I Pi + "^^ pi + "^^ 9i + "^^ ^2 + '^ \ 

+ p + 2)iVi?(p) ^^1 +^2 + 'Zi + 'Zs) + (^ + 2)iV ^ 

NB^ip) 
where we introduced a second bare coefficient 



y(°) = V / -^ '^J^ ^ = Y + 0{e) (B 46) 



The formula (B.45) is now ready for renormalization and for the renornialized four-point 
function we get (2.48) of section 2. 

It remains to calculate the finite part Y . We proceed similarly to the evaluation of the 
coefficient oi X^°' above but here the calculation is more difficult because of the momentum 
dependence. We start by defining 






(B.47) 



{27t)D (g2+j„2)2[^2 + (g+p)2] 

+ ^2iSal3P-yP6 + ^a-yP/SPS + ^aSPjBP-i + PaP 13^-^5 + PaP^y^/SS + PaPS^/S-y) 
+ ^3PaPl3P-fPS- 

What we need is 

Y / TT^T-^ ^?v?r^9 7 7^=3I?ni+6fi2/ + r?3VPu (B.48) 

^J {2tt)^ {q^ + m?Y[m? + {q + pY] •^Z_^^/^ ^ ' 

and we expect that 

n = 3DJ]i + 6J72p'--— ^-— =a; + 0(e), J^s = ^3 + O (e) . (B.49) 

{U + 2) Jo 

Anticipating the finiteness of tJ and co^ we can write 

Y = lJ + uj3p'^. (B.50) 
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For this calculation we have to introduce, in addition to (B.47), the following hierarchy 
of tensor integrals (r = 0, 1): 



+ NrPaPlSP-r, 

f d^g ^^^^^ (B.51) 

/ 



ZrPa- 



{2tt)D (g2 + m2)(i+'^) [m2 + (g + ;p)2] 

By making contraction of the indices with Kronecker deltas and the momentum vector p^ 
we find the following simple relations among the tensor integrals: 

{D + 2)^1 + p'^n2 = Ho-m'^Hi, (L> + 4)f^2 +/f^3 = Mo-m'^Mi, 

2!^i + 2p^n2 = -Kq - p^Ki, 6JI2 + 2/!^3 = -No - p^Ni, 

{D + 2)Ko+p^No = -I//0 - m%, {D + 2)Ki +p^Ni = Co - m^^i, 

2/K0 = V^o, 4Ko + 2p^No = -I//0 - /Mo, 

2p^Ki = (I//0 - m25(0))/Z) -Ho- p^H^, AKi + 2p2Ar, = -Mo - p^M^, 
DHo+p^Mo = I//0 - m^B{p), DHi+p^Mi = B{p) + {rr? /2){dB{p)/dm^), 

2Ho + 2/ Mo = I//0 - /Co, 2Fi + 2p2Mi = -^o - p'^i, 

2p^eo = -p^Bip), 2/^1 = 5(0) - B{p) + (pV2)(95(p)/9m2). 

(B.52) 
We have solved this overdetermined system of linear relations using Mathematica. The 
results for ilj are too bulky to be reproduced here. In the limit D — )• 2 the formulas 
simplify enormously and we find 

^ = -3™=.(rt + ^ U + ^1 ?^ + i f 5 + 4 V (B.53) 



ij,'M = 2mH(j.) -{4 + ™V + ^"l ?^ - T^ (5 + 4 1 ^ (B^54) 



4 \ 4 J dm? 32tt \ m? 

2 ^"^P ^ 4 j dm? SvrV m2 
C Lattice integrals 

First note that the continuum limit of the lattice integral -Biatt(p) in eq. (4.15) is given by 
h{p) defined in eq. (2.52) which can be expressed analytically as 



1 / yp^ + 4m2 + 



'p 



2 



b{p) = b{p, m) = „ , ,, , ^ - , In \" . . ^P ^ , (C.l 



2tx\Jp^{jP' + 4m2) \ y/p2 _|_ 4j^2 

where the explicit reference to the second argument m will be used later in this appendix. 
6(p) is analytic in i? with a cut from —00 to —4rr?. Also note that 6(p) 7^ for all ]? and 
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the properties: 



6(0) 



1 



Further we have 



d 



dm? 

d 

dPf, 



b{p) 
h{p) 



AiTm? 
2 



p^ + Am? 



Kp) + 



1 



Airm? 
(/ + 2m^)h{p) 



1 

2^ 



(C.2) 
(C.3) 

(C.4) 
(C.5) 



In order to work out the leading lattice artifacts of -Biatt {p) it is convenient to consider 
the equivalent representation: ^^ 

2 



Aatt(p) = a^ dtidt2 e-(*i+*^)('^'™o+4) Jl Jo (^2^(ti + t^f - ht2a^p] 

/■oo />! 2 

= a^ / dt / dxte-*("''"o+4) TT j^ (2tS^(x,p)) , 
Jo Jo ^^1 



(C.6) 



where 



SfM{x,p) = Y 1 - 2;(1 - x)a2p2 ^ 
and /^ denote the modified Bessel functions. For p = we have simply: 

POD 

Aatt(O) = a^ / dtte-'^^'-o* [e-2*/o(2t)]' . 
Jo 

Noting for large t: 



e-^'lri2t) 



'A-Kt 



1 _ 4r2 

i + i^ + e.(.-) 



(C.7) 
(C.8) 
(C.9) 



and splitting the range of the t-integration in (C.8) into parts [0, 1] and [l,oo] we obtain: 



Aatt(O) 



AirmQ 



+ a^ 



ozvr 



+ 0{a^mlln{a^ml)) 



with 



ci 






(fc) 



fc=i 



f> = y"d(([e-''/„(2()]=, 
f'^fd.{ne-%(2.,]'-^(l.l)}, 



.(3) 



.(4) 



1 



327r 



dtt-^e"*. 



(CIO) 

(C.ll) 

(C.12) 
(C.13) 
(C.14) 
(C.15) 



"^before performing the g^ integrations we shift the variables g,i -^ q^^ + a^ with sin(aQ^) = 



-i2 sin(app)r^ ^ , cos(aQfi) = [ti + t2 cos(ap^)] T^^ ,T^, = i/(fi +f2)^ - tit2a'^p- 
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ci can in fact be determined analytically as 



ci 



1 
16^ 



-ln(2)-l 



which can easily be checked numerically. 
Proceeding as for the case p = above: 



-Biatt(p) = a^cY +a 



/ dt / dxte-*('^'"o+4) TT iQ(2tS^{x,p)) + O (a^ 



1 TOO 



+ o^-— / dt / dx' 



-t[a2m2+4-2E'=i5.(x,p) 



1 + ^ 



+ O (a^) 



c(^) + cf ) 



+ 



47r 







X <! 1 + ^ + ^a;(l - x)aV - y^ 



(l-x)2aV|+0(a') 



6(p, mo) + a^ 



ci + ^/3k 



k=l 



+ O {a') , 



with 



/3i 



'dx. '^'-'^P" 



Wit Jq [m^ + x(l — x)^^] 

x^(l — x)^p^ 



1 1 2,. ^ 



1 /-l 

/32 = / dX ■ r, 

iDTT 7o [m2 + x(l -X)p2]^ 



p-^ 



1 5 
2inn?b{p) - m^- — Mp) 



Ps 



4(p2)2 

— / dx In Ta m + x(l — x)a p 
27r Jo 



2„21 



327r 
1 
16 



{p^ + 4m^)b{p) H ln(a m ) 

2tt tt 



(C.16) 



(C.17) 



(C.18) 



(C.19) 
(C.20) 



(C.21) 
(C.22) 
(C.23) 

(C.24) 
(C.25) 
(C.26) 



Putting the above results together and expanding rriQ and p in m and p respectively 
(i.e. in a) we obtain (4.16) with 



hi{p) 



m 



4 d 



12 97^2 



^(^) - E ft ^^(P) + ci + /3i + /32 + /33 



(C.27) 



Using eqs. (C.4),(C.5) and inserting the resulting bi{p) into eq. (4.23) one indeed obtains 
eq. (4.24). 
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For the tadpole (which we need below) one gets 



J(mo) = 



w/a j2 



1 



-TT/a (2vr)2 ml + q 



Air 



ln(a itiq) + C2 + a ttt-q 



327r ^ °^ 327r 



where 



., = |_,n(2)=8c. + i-. 



(C.28) 
+ O (a^) , (C.29) 

(C.30) 



Next we consider the integrals occurring in the 4-point function of the mixed action. 



First 



VW 



TT/a ^2q 



Qt^(.P + Q), 



n/a (2vr)2 



(mg + f) 



ml + {p + q) 



2cos (^) 1^ dt J' dxte-'^-"^o+^^Io{2tSj,{x,p)) 
h{2tS^{x,p)) 



(C.31) 



/o ^0 

/i(2t5^(x,p)) 



where 1 = 2,2 = 1. The expansion for /i^ (0) to order a? is immediately obtained from 
eq. (C.29): 



^(0) 



d 

J{mo) + ml-—^J{mo) 



1 

8tt 



2 Svr 



"l^2 2\,'-2 -^,2 2 

— ln(a mQ) + — - — + amo 



^ln(a^mg)-^-ci 



(C.32) 



+ 0(a^) . 



Actually for our purposes we only need /i^(p) to leading order 

h^ip) = 2h^^\p)-plh(;^Hp) + 0{a^) , 

/"OO /"l 

h^^){p)= dt / dxte-*('*''"o+4)jQ(2tS-(x,p)) 

Jo JO 

x[/o(2t5^(x,p))-/i(2tS^(x,p))] , 



/^J.^)(P) 



OO 



a- / dt / dxtx(l - a;)e-*(" ™o+4)/q (2t5^(x,p)) /i (2t5/,(x,p)) . 
'0 Jo 



From these representations we can deduce the leading terms 

where the (3^ are defined in eqs. (C.21)-(C.25). Adding the terms we get 



^m(p) 



2 

^^ + ^-W + ^^'Mp) + ^b{p)Q + O (a2) 
2 ovr 4 p^ ^ ' 



(C.33) 

(C.34) 

(C.35) 

(C.36) 
(C.37) 

(C.38) 
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with (5, z defined in eqs. (2. 53), (1.19) respectively. 
Finally for 



^ViviV) 



--/'^ ^^""^^ {ml + g2) \ml + (f+ q) 



2 



(C.39) 



appearing in the definition of /\^u{p) (cf. [8], app. B) and which we need here only for 
p = 0, mo = 0, we have 
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